Dynamical instability of the XY spiral state of ferromagnetic condensates 
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We calculate the spectrum of collective excitations of the XY spiral state prepared adiabatically 
or suddenly from a uniform ferromagnetic F = 1 condensate. For spiral wavevectors past a criti- 
cal value, spin wave excitation energies become imaginary indicating a dynamical instability. We 
construct phase diagrams as functions of spiral wavevector and quadratic Zeeman energy. 



Spinor condensates of ultracold atoms [H, [H, d, 0, @, 
0, H are the latest addition to the class of many-body 
systems with multicomponent order parameters. One of 
the most intriguing manifestations of the high symmetry 
of such systems is the possibility of a large variety of spin 
textures and topological defects. While earlier studies of 
liquid crystal nematics and superfluid 3 He demonstrated 
the existence of spin textures, experiments with spinor 
condensates provide a unique opportunity to investigate 
their non-equilibrium quantum dynamics. Understand- 
ing dynamical properties of spin textures will provide 
valuable insight into many open problems of quantum 
magnets and spinful superfluids, from analysis of the 
Kibble-Zurek mechanism of nucleating topological de- 
fects when crossing a quantum phase transition [Io| to 
finding the fundamental limits of spinor BEC magne- 
tometers 
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In this paper we investigate theoretically the stabil- 
ity of the spin spirals in ferromagnetic S=l condensates 
(see Fig. 1.) Such states represent the simplest type of 
spin structures and can be prepared experimentally by 
applying a gradient of the magnetic field in the direction 
perpendicular to the magnetization axis [12j |. Our main 
result is the prediction of dynamical instabilities for spiral 
states, which we summarize Fig. 2. We find that rotation 
of the magnetization vector from the XY plane to the z 




FIG. 2: (color online) Collective mode phase diagrams for 
fz = against spiral wavevector k and quadratic Zeeman 
energy q in the adiabatic and sudden (inset) limits. DLi (DI[) 
indicates a dynamical instability with one branch of unstable 
modes beginning at k = (k > 0). DI2 indicates a dynamical 
instability with two distinct branches of unstable modes. 



direction drives the instability at small quadratic Zeeman 
energy while rotations within the XY plane are respon- 
sible at large magnetic fields. Surprisingly we observe 
that the unstable modes are characterized by wavevec- 
tors that can be considerably larger than the wavevector 
of the initial spiral state. Other instabilities in spinor 
condensates discussed previously include Castaing insta- 
bilitie s 11311 in incoherent non-condensed two component 
87 Rb [14| and the modulational instability of a uniform 
spinor condensate ljl 1|| 17 1. 

Our starting point is the microscopic Hamiltonian 



H = * f 
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FIG. 1: (color online) From left to right: magnetization vec- 
tor in the XY spiral state for fully polarized, partially polar- 
ized, f z 7^ 0, and after fragmentation. 



where \f a with a = x, y, z are annhilation operators for 
F = 1 bosons with mass m and (F a )^ are angular mo- 
mentum operators. We use a matrix notation with sup- 
pressed indices where *, T, and f denote the complex 
conjugate, transpose, and the conjugate transpose, re- 
spectively. For example, ^ (Nfl) is a column (row) vector 
while F z is a matrix 
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Interaction strengths are given by go = Airffao/m, 
g s — ATrh 2 (ao — 02) /3m Q in terms of the s-wave scatter- 
ing lengths a p for two atoms colliding with total angular 
momentum F and : : denotes normal ordering. For 87 Rb, 
a = 101. 8as and 02 = lOOAas where as is the Bohr 
radius [l8| giving positive g s and ferromagnetic interac- 
tions. 

This Hamiltonian has a U(l) (8 SO (2) symmetry of 
global phase rotations and spin rotations about the z 
axis. The chemical potential fi and linear Zccman energy 
p are Lagrange multipliers controlling the corresponding 
conserved quantities 



n, 



(*tF,*) - nf„ 



(2) 



where n is the total particle density and f z is the z com- 
ponent of the magnetization per particle. 

Due to conservation of F z , static magnetic fields enter 
through the quadratic Zeeman energy q instead of the 
linear Zeeman energy p. Moreover, q can be further ma- 
nipulated through the AC Stark shifts. From here on, we 
take representative values q = 70 Hz G~ 2 B 2 where B 
is the magnetic field and n = 2.2 x 10 14 cm" 3 We 



neglect here magnetic dipole interactions 

The XY spiral state is prepared from an initial cigar 
shaped condensate with uniform XY magnetization by 
applying a magnetic field gradient along the axial or z 
axis. After switching off the gradient, imaging of the 
transverse magnetization can be used to study the sta- 
bility of the XY spiral state. As illustrated in Fig. [TJ 
the transverse magnetization winds along the z axis. The 
transverse magnetization is fully polarized for q = f z = 
and is suppressed due to population of the m z = com- 
ponent of \& when q 7^ or f z 7^ 0. In the presence of an 
instability, the system fragments into domains carrying 
different magnetization vectors. 

Analyzing the generation of the XY spiral state re- 
quires understanding of the complicated non-equilibrium 
dynamics. We focus on studying the resulting non- 
equilibrium stationary state which we assume is well- 
described as a coherent condensate. Such states are given 
by mean- field solutions of the Gross-Pitaevskii (GP) 
equations implied by Eq. [T] which carry XY spiral order. 
Compared to stable ground states, these non-equilibrium 
stationary states are in general metastable and decay via 
linear and non-linear processes. We consider the linear 
stability of such states with respect to small fluctuations 
by analyzing the spectrum of collective modes. The dis- 
tinction between stable and metastable stationary states 
also arises for spinless bosons in a moving optical lat- 
tice 2l|, |22j and in optics in the context of the four-wave 



mixing instabilities (like the superradiance instability), 
in which a mean-field treatment may yield a stationary 
situation only because it neglects the spontaneous scat- 
tering into initially unoccupied modes of the system Q • 
To perform the stability analysis, it will be useful to 
consider a frame comoving with the XY spiral order. We 



thus introduce the substitution 



* — * exp (inzF z ) * 



(3) 



where k is the spiral wavevector. The comoving frame 
Hamiltonian is then given by Eq. [T]with the subsitution 



p^p + 



IK. 

— v„ 

m 



<l 



2m < 4 > 



We use this spin-dependent effective Hamiltonian to 
study the stability of the non-equilibrium XY spiral state. 
However, we note that it may be possible to explicitly en- 
gineer a physical Hamiltonian of this form through con- 
tinuous Raman excitation similar to that described in 
Ref. 0. 

In the adiabatic limit which we describe first, the com- 
ponents of '3/ are able to adjust their populations in or- 
der to accomodate the XY spiral order. The interaction 
terms of Eq. [T] describe spin flip processes that mix the 
components of ^ but still conserve the overall magneti- 
zation. The components of St in the resulting XY spiral 
state then describes a compromise between the kinetic 
energy cost of the winding spiral and gain in interaction 
energy. 

We begin by looking for mean-field solutions of the 
GP equations in the comoving frame of the form \I f = 
y/nffre^ 1 where we use the following parametrization 



e"V 1 -/ Z coth (2x) 



(5) 



which automatically gives the correct conserved quanti- 
ties of Eq. [2] by construction. 

The parameter rj describes a global phase that spon- 
taneously breaks U{1) phase rotation symmetry of the 
Hamiltonian. Similarly, (j) gives the orientation of the 
magnetization vector and breaks SO(2) spin rotation 
symmetry. Here rjj_ gives a relative phase between the 
z and transverse components of \&. Solutions of the GP 
equations for r]± distinguish between g s > ferromag- 
netic and g s < antiferromagnetic interactions where 
r]± = and rj± — Tt/2, respectively. Recall we focus on 
the g s > case. Finally, \ controls the relative magni- 
tude between the z and transverse components of '. Af- 
ter introducing the dimcnsionlcss quantities r = tanh(x) 
and Q — q/2g s n we find the GP equations give the con- 
dition 



Qt 3 



[l-Q)r = f z 



(6) 



As in Ref. [25| , we find three classes of mean- field solu- 
tions: polar, ferromagnet, and XY spiral state. The polar 
state occurs for f z — and Q > 1 while the ferromagnet 
occurs for f z = ±1. Both of these states occur only on 
isolated lines in the mean-field phase diagram and do not 
support XY spiral order. 
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FIG. 3: (color online) Representative collective mode dis- 
persions for f z = and q = 0.2 Hz in the adiabatic limit 
illustrating one (two) branches of unstable modes, c (s) de- 
notes the charge (spin) mode. 



We now analyze the stability of the XY spiral state by 
studying the spectrum of collective fluctuations 6<f> about 
the mean-field solution <&. We take '8/ = ^/n(<&+5&)e lu ' t . 
and focus first on the f z = case. Using the standard 
Bogoliubov analysis, we find the excitation energies ujk 
satisfies the eigenvalue equation 



det 



M k - uo k N 

-N* -M*_ k - u k 







(7) 



where and N are given by 



M k = 



P 



2m \ m ) \ Ira 



F 



(8) 



Recall we use a matrix notation so that <I>^<1> (3><I>f) is 
a scalar (matrix). We consider the one-dimensional case 
relevant for cigar shaped condensates where we can take 

The XY spiral state spontaneously breaks U(l) <g> 
50(2) symmetry of global phase and spin rotations and 
we find both a gapless charge and spin mode (as required 
by the Goldstone theorem) with linear dispersions. How- 
ever, the spin mode can develop imaginary frequencies 
which indicate the presence of a dynamical instability. 

We find there are several distinct types of behavior for 
the spin mode exhibiting a dynamical instability. The 
first possibility is a branch of unstable modes starting at 
k = which we denote as DI\ as illustrated in Fig. [3] 
The second is a branch of unstable modes starting modes 
starting at k > which we denote as DI[ . The third pos- 
sibility is two distinct branches of unstable modes start- 
ing at k = and k > as illustrated in the inset of Fig. 

El 

We construct the phase diagrams illustrated in Fig. [2] 
by characterizing the behavior of the spin mode as a func- 



tion of the spiral wavevector k and quadratic Zeeman en- 
ergy q. We first consider the adiabatic limit characterized 
by an interpolation between long-wavelength instabilities 
in the limit of large and small q. Both instabilities can be 
throught of as unwinding of the spiral order through ro- 
tations of the magnetization vector, but they arise from 
qualitatively distinct physics. 

When q is zero, the system is rotationally symmetric. 
In this limit, the XY spiral state is potentially unstable 
towards unwinding through arbitrary S*0(3) rotations of 
the magnetization vector from the XY plane to the z 
axis. However, small but finite q provides a potential 
energy barrier to such a process. When the kinetic energy 
stored in the non-uniform winding is sufficiently large, 
small fluctuations corresponding to such rotations can 
overcome this energetic barrier and grow exponentially 
giving rise to a dynamical instability. In particular, the 
instability in this regime is due to fluctuations in the 
direction of the magnetization vector. 

In contrast, large q explictly breaks rotational sym- 
metry. The magnetization vector is essentially confined 
to the XY plane and the unwinding of the XY spiral 
state can then only proceed via SO{2) rotations within 
that plane. Such rotations proliferate near the quantum 
phase transition to the polar state when fluctuations in 
the magnitude of the magnetization vector are large. 

This large q instability can be mapped to the insta- 
bility of current carrying states for spinless bosons [2t| . 
Here the effective SO(2) magnetization order parameter 
maps to the U{1) order parameter of spinless bosons. In 
addition, the critical fluctuations near the transition to 
the polar state driving the instability map to those of 
bosons near the Mott transition. 

From the physical arguments above, we expect the 
XY spiral state to be stable for wavevectors less than 
k 2 /2m ~ q for small q when kinetic energy is insuffi- 
cient to overcome the potential energy barrier. Similarly, 
the XY spiral should be stable for wavevectors less than 
k 2 /2m ~ (q — q c ) near the quantum phase transition to 
the polar state at q r . The boundaries in Fig. [2] can be 
obtained explicitly [27j as 



2m 



< 



2g s n - q 



2^ 

90 



q> 



2m 



(9) 



which gives k 2 /2m < q and n 2 /2m = (q c — q)/3 for for 
small and large q, respectively, in agreement with the 
physical arguments. 

Also notice in Fig. [2] an isolated line on which the 
XY spiral state is stable at small B and intermediate n. 
In fact, the dynamical instability in the DI[ surrounding 
this line is weak in the sense that the imaginary part of oju 
is comparatively small. The energetic arguments given 
earlier for the small B limit seem to suggest increasing k 
makes the XY spiral state more unstable. However, as n 
increases further, the populations in the components of 
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FIG. 4: (color online) Collective mode phase diagrams for 
f z = 0.05 (left) and f z — 0.5 (right) against spiral wavevec- 
tor k and quadratic Zeeman energy q in the adiabatic (top) 
and sudden (bottom) limits. Dl (LI) indicates a dynamical 
(Landau) instability. DI\ , DI[ , and DI2 are described in Fig. 
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\& change appreciably. 

In this limit, the XY spiral state can thought of as 
having a significant polar component. It has been shown 
previously that the polar state has a dynamical instabil- 
ity with a characteristic wavevector k ~ ^/2mg s n [l5j |. 
When the spiral wavevector is on the order of this char- 
acteristic wavevector, the XY spiral state becomes less 
unstable. Physically the spiral order tends to suppress 
the instability of the polar component. 

After analyzing the f z = case in detail, we now briefly 
discuss the f z ^ case. The collective mode phase dia- 
grams for the adiabatic limit are shown in the top of Fig. 
Hwith f z = 0.05 to the left and f z = 0.5 to the right. 
The small f z phase diagrams are qualitatively similar to 
the f z = case. However, there is no polar state which 
only occurs for f z — but there is an additional region 
which exhibits a spin mode exhibiting a dispersion with 
negative frequencies, corresponding to a Landau insta- 
bility. For large f z , the phase diagrams no longer exhibit 
a characteristic peak for the stable region as a function 
of q. 

So far the results have been focused on the adiabatic 
limit where the components of can adjust due to mag- 
netization spin flip processes. In practice, the prepara- 
tion of the XY spiral state can also occur on a timescale 
shorter than that of spin flips. We thus briefly comment 
on qualitatively similar results in the sudden limit where 
the populations of each component cannot change from 
their initial values. To take this effect into account, we 
consider mean-field solutions of the GP equation of the 
form = ^/n^e lUat with u> x = u> y ^ uj z and 3? is as in 
Eq. [5] Notice the components of \& evolve at different 
frequencies which allows for solutions with the necessary 



populations for each component. 

We then perform the same analysis of the collective 
modes as in the adiabatic limit. This gives for f z = in 
the sudden limit the phase diagrams in the inset of Fig. 
[2] As in the adiabatic limit, the sudden limit is charac- 
terized by an interpolation between long-wavelength in- 
stabilities in the small and large B limits arising from the 
same physical origins. The region where the XY spiral 
state is stable is given by 



< 



2g s n 



2m ~ 2+2^ ' 

90 




(10) 



which gives k 2 /2m = q/2 and k 2 /2m = (q c — q)/2 for the 
small q and large q limits, respectively. Up to coefficients, 
this is of the same form as the adiabatic limit. The f z = 
0.05 and f z — 0.5 phase diagrams for the sudden limit 
are shown in the bottom of Fig. [4] and exhibit the same 
structure as the adiabatic limit. 

In this paper we have focused on the one-dimensional 
limit relevant for cigar shaped condensates. However, 
the formalism we used can be readily adapted for the 
three-dimensional case. In particular, one just needs to 
take k z = k cos 9 where 6 is the angle between the mode 
wavevector and spiral wavevector in Eq. [8j 

In summary, we studied a possible mechanism for the 
instability of the XY spiral state. Focusing on the limits 
where the XY spiral is prepared adiabatically or sud- 
denly, we demonstrated that when the spiral wavevector 
exceeds a critical value spin wave energies become imag- 
inary. This indicates the presence of a dynamical insta- 
bility and exponential growth of fluctuations. We traced 
the physical origin of these instabilities to unwinding of 
the magnetization vector through rotations from the XY 
plane to the z axis for small quadratic Zeeman energy q 
and within the XY plane for large q. 
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uate Research Fellowship, Harvard-MIT CUA, AFOSR, 
MURI, and the NSF grant DMR-0132874. When this 
work was being completed we learned of the work by A. 
Lamacraft addressing related issues [28| . 
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